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Abstract 

A search for stable and long-lived massive particles of electric charge |Q/e| = 1 or 
2/3, pair-produced in e^e~ collisions at centre-of-mass energies from 130 to 183 GeV, is 
reported by the OPAL collaboration at LEP. No evidence for production of these parti- 
cles was observed in a mass range between 45 and 89.5 GeV. Model-independent upper 
limits on the production cross-section between 0.05 and 0.19 pb have been derived for 
scalar and spin-1/2 particles with charge ±1. Within the framework of the minimal su- 
persymmetric model (MSSM), this implies a lower limit of 82.5 (83.5) GeV on the mass of 
long-lived right- (left-)handed scalar muons and scalar taus. Long-lived charged leptons 
and charginos are excluded for masses below 89.5 GeV. For particles with charge ±2/3 
the upper limits on the production cross-section vary between 0.05 and 0.2 pb. All limits, 
on masses and on cross-sections, are valid at the 95% confidence level for particles with 
lifetimes longer than 10~^ s. 
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1 Introduction 



Most searches for new particles predicted by models beyond the Standard Model (SM) assume 
that these particles decay promptly at the primary interaction vertex due to their very short 
lifetimes. These searches would not be sensitive to long-lived heavy particles which do not 
decay in the detectors. However several models predict such long-lived particles. For example, 
in the minimal supersymmetric model (MSSM) if the mass difference between the chargino 
and the lightest neutralino is smaller than a few hundred MeV, the lightest chargino would have 
a lifetime sufficiently long to result in decays predominantly outside the detector. In gauge- 
mediated supersymmetry, if the slepton is the next-to-lightest supersymmetric particle, it could 
decay with a long hfetime to a lepton plus a goldstino, if the SUSY-breaking energy scale is 
sufficiently high 0. R-parity violating SUSY models also allow for long-lived heavy charged 
particles. If the lightest supersymmetric particle is a slepton and the R-parity violating coupling 
is small (A < 10~®), then the slepton decay length would be larger than a few meters, the typical 
size of tracking detectors at LEP. If a fourth-generation heavy lepton (mL > mzl'2) exists g], 
the charged heavy lepton would be stable if it is lighter than its neutral isodoublet partner. 
Some models beyond the SM would also predict the existence of particles with fractional electric 
charge. As an example, leptoquarks having a small coupling Al,r could be long-lived and 
possess fractional charge. Another example could be long-lived hadronic states with fractional 
charge predicted by some modified QCD models 0. 

Previous searches for long-lived charged particles have been performed by the LEP collabora- 
tions with data taken at the TP resonance |]^. The OPAL Limits on the production cross-section 
have been set at approximately 1.5 pb at 95% CL for masses between 34 and 44 GeV. DELPHI 
and ALEPH |^ have also analysed the data collected at centre-of-mass energies up to 172 GeV, 
setting limits of 0.2-0.4 pb for masses between 45 and 86 GeV. 

This paper describes a search for long-lived particles X^, with mx > with charge 

\Q/e\ = 1 or 2/3, pair-produced in the reaction e"'"e~ — > X"''X~(7). The data were collected 
by the OPAL detector during 1995-1997, at centre-of-mass energies of 130-136, 161, 172 and 
183 GeV for a total integrated luminosity of 89.5 pb~^. Such particles distinguish themselves 
by their anomalous, high or low ionization energy loss, dE/dx, in the tracking detector gas. 
This search is therefore primarily based on the precise dE/dx measurement provided by the 
OPAL jet chamber. However, for particles of charge one, there is a large mass region where the 
measured dE/dx cannot distinguish them from ordinary particles. A complementary search, 
based on the two-body kinematics, is used to cover this mass region. No search was made 
for particles with \Q/e\ = 1/3 because their low ionization energy loss is too close to the jet 
chamber dE/dx measurement threshold. The results obtained are valid for particles with a 
lifetime longer than 10~^ s. 

2 The OPAL Detector 

A complete description of the OPAL detector can be found in Ref. and only a brief overview 
is given here. The central detector consists of a system of tracking chambers, providing charged 
particle reconstruction over 96% of the full solid angleQ inside a 0.435 T uniform magnetic field 

^The OPAL right-handed coordinate system is defined such that the z axis is in the direction of the electron 
beam, the x axis is horizontal and points towards the centre of the LEP ring, and 9 and (j) are the polar and 
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parallel to the beam axis. It consists of a two-layer silicon microstrip vertex detector, a high- 
precision drift chamber, a large-volume jet chamber and a set of 2;-chambers measuring the 
track coordinates along the beam direction. 

The jet chamber (CJ) is the most important detector for this analysis. It is divided into 24 
azimuthal sectors, each equipped with 159 sense wires. Up to 159 coordinate and dE/dx mea- 
surements per track are thus possible, with a precision of ar,j, ~ 135 /zm and ~ 6 cm. When 
a track is matched with 2;-chamber hits the uncertainty on its Zq coordinate is ~ 1 mm. The 
jet and 2;-chambers, located inside the magnetic coil, provide a track momentum measurement 



with a resolution of a.p/p ~ y (0.02)^ + (0.0015 ■ pt^ for tracks with the full number of hits {pt, 
in GeV, is the momentum transverse to the beam direction) and a resolution on the ionization 
energy loss measurement of approximately 2.8% for fi~^fi~ events with a large number of usable 
hits for dE/dx 



A lead-glass electromagnetic calorimeter (ECAL) located outside the magnet coil covers 
the full azimuthal range with excellent hermeticity in the polar angle range of | cos^^| < 0.984. 
The magnet return yoke is instrumented for hadron calorimetry (HCAL) covering the region 
I cos^l < 0.99 and is followed by four layers of muon chambers. Electromagnetic calorimeters 
close to the beam axis complete the geometrical acceptance down to 24 mrad on both sides 
of the interaction point. These include the forward detectors (FD), which are lead-scintillator 
sandwich calorimeters and, at small angles, the silicon tungsten calorimeters (SW) |Tl|. The 
gap between the endcap ECAL and the FD is instrumented with an additional lead-scintillator 
electromagnetic calorimeter, called the gamma-catcher (GC). 

The ionization energy loss dE/dx produced by a charged particle is a function of (3^ = p/m 
and of the electric charge Q |]lOl- In Figure |l|, the distribution of dE/dx as a function of the 
apparent momentum, p/Q, is shown. Standard particles of charge ibl (e, /x, tt, p, K) with 
high momentum (p > O.l-y/s GeV) have dE/dx between 9 and 11 keV/cm. Massive particles 
with charge ±1 are expected to yield dE/dx > 11 keV/cm for high-mass values, mx > 0.36^^ 
(e.g., at y/s = 183 GeV, mx = 65 GeV, one has p = QA GeV and dE/dx = 11 keV/cm), or 
dE/dx < 9 keV/cm for low-mass values, mx < 0.27^^5 (e.g., at i/i = 183 GeV, mx = 50 GeV, 
one has p = 77 GeV and dE/dx = 8 keV/cm). The dE/dx measurement therefore provides 
a good tool for particle identification in this high- and low-mass regions. However, in the 
intermediate mass region, the dE/dx measurement does not distinguish the signal from the 
background. In this region an analysis based on kinematic properties of pair-produced massive 
particles is used. 



3 Monte Carlo Simulation 

This section describes the Monte Carlo simulation of the signal and the background samples. 
All generated events have been processed through the full simulation of the OPAL detector |T2]; 
the same event analysis chain has been applied to the simulated events and to the data. 

To generate the signal process e"'"e~ X+X^(7) different generators have been used. Signal 
events of the type e"'"e~ — > (7) (£=*= being a scalar lepton) have been generated at four 
different energies (a/s = 133,161,172,183 GeV) using SUSYGEN |jl3[. The generated scalar 



leptons are not allowed to decay, therefore simulating the signal from heavy charged stable 



azimuthal angles, defined relative to the +z- and +x-axes, respectively. The radial coordinate is denoted by 
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scalar particles. Similarly, events of the type e^e" L+L^(7) and e^e" QQ(7), where 
are stable heavy leptons and Q, Q are stable heavy deconfined quarks with charge ±2/3, have 
been generated at the same energies, using the generator EXOTIC |T^. All signal samples have 
been generated with mx ranging from 45 GeV to 90 GeV. Each sample contains 1000 events. 
For the purpose of detector simulation and particle interactions, and L''^ particles have been 
treated as heavy muons, while Q, Q as heavy stable hadrons with charge ±2/3. 

The background has been estimated using simulations of all Standard Model processes 
(lepton-pair and multihadronic processes, four-fermion processes, two-photon processes). The 
Monte Carlo samples generated at -^i = 171 and 184 GeV are briefly described below. Small 
differences in the centre-of-mass energies between data and the simulated background Monte 
Carlo samples have a negligible effect on the analysis. A detailed description of the Monte 
Carlo samples generated at ^/s = 133 GeV and = 161 GeV can be found in |15|, |1^. All 
background samples have an equivalent luminosity of at least ten times the data collected at 
each energy. 

The contribution to the background from two-fermion final states has been estimated using 
BHWIDE [0 for the e+e-(7) final states and KORALZ g for the /i+/i-(7) and the r+r-(7) 



states. Multihadronic events, qq(7), have been simulated using PYTHIA [ p!9[| . 

For the two-photon background, the PYTHIA [0, PHOJET ^ and HERWIG ^ Monte 
Carlo generators have been used for e"'"e~qq final states and the Vermaseren generator for 
all e'^e~i~^i~ final states. All other four-fermion final states have been simulated with grc4f ||23 
which takes into account interferences between all four-fermion diagrams 



4 Data Analysis 

Pair-produced stable or long-lived massive charged particles would manifest themselves as 
events with two approximately back-to-back charged tracks. Particles of charge ±1 are as- 
sumed to not interact strongly, hence to not produce hadronic showers. Since they are massive, 
they do not produce electromagnetic showers either. From these considerations, these events 
should be very similar to fi~^^~ events, the only difference being the higher mass of the particles. 

No assumption is made for the interaction properties of long-lived particles with charge 
±2/3; hence no calorimetric signature is used in the search for such particles, which is described 



in subsection 4.3 



The events collected by OPAL at y/s =130-183 GeV have been required to pass a preselection 
which rejects events incompatible with the signal topology. The following criteria are applied: 

PI Events are rejected if the total multiplicity of tracks in the central detector and clusters 
in the ECAL is greater than 18. Cosmic ray events are rejected as well as Bhabha 
scattering events . 



P2 Events are required to contain exactly two tracks in the central detector satisfying basic 
quality criteria^ and having a momentum p > 0.1^/s, a momentum transverse to the 



^The distance between the beam axis and the track at the point of closest approach (PCA) must be less 
than 1 cm; the z-coordinate of the PCA must be less than 40 cm; the innermost hit of the track measured by 
the jet chamber must be closer than 75 cm to the beam axis. 
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beam axis pt > 0. 025^/5, polar angle satisfying | cos 6*1 < 0.97 and at least 20 CJ hits 
usable for dE/dx measurement. The two selected tracks are required to have opposite 
electric charge. 

P3 To reduce background from two-photon interactions, the acoplanarity angleQ between the 
two tracks is required to be 0acop < 20° and the total visible energy^ of the event is 
required to be E^is > 0.2y/s. 

P4 To reduce background from events with initial state radiation, events containing an iso- 
lated ECAL cluster with an energy greater than 5 GeV are rejected. Isolation is defined 
as an angular separation of more than 15° from the closest charged track. 

P5 It is required that ^ + ^ < 0.2, where £'1^2 denotes the energies of the ECAL clusters 
associated to the two selected tracks, to further reduce the contribution from Bhabha 
scattering events. Moreover, in a cone of 10° half-opening angle around each of the two 
selected tracks, no other tracks with p > 0.5 GeV and no unassociated clusters with 
E > 3 GeV should be found. 

After these preselection criteria, the background is dominated by e"'"e~ — > events, 
with a small contribution from e^e~ t^t~ and two-photon e+e~/i"'"/i~ events. The effect of 
these preselection cuts on the samples at \/s = 183 GeV is shown in Table 0. 

4.1 Search for particles with charge ±1 

Two complementary methods are adopted, depending on the mass rrix of the signal particle. In 
the high-mass (mx > 0.36y^) and low-mass (mx < 0.27y^) regions, the dE/dx measurement 
is used to distinguish the signal from the background. The preselected events are retained if 
they satisfy the following requirements on dE/dx: 

Al Both high- momentum tracks must have either dE/dx > 11 keV/cm or dE/dx < 9 
keV/cm. 

A2 The probability that either of the two dE/dx measurements originates from one of the 
standard particles (e, /i, vr, p, K) must be less than 10%. 

After this selection (A1-A2), no candidate survives in any of the data sets. The total 
background is estimated to be 0.03 events at ^/s = 130 — 136 GeV and less than 10~^ events 
at any other energy (see Tables |l|, |^). 

A complementary analysis is used for masses in the intermediate mass range (0.27 < 
"^x/a/s < 0.36), where the dE/dx measurement does not provide adequate separation be- 
tween the signal and the di-lepton background. In events of the type e^e~ — >■ X+X~, since 
mx > mz/2 is assumed, initial state radiation is suppressed. Hence, each particle should have 
an energy close to the beam energy and nix can be estimated for each track as mx = \J s/A — p"^. 
However, for SM events with two tracks and missing energy (e.g., events with neutrinos, with 

•^The acoplanarity angle, 0acop, is defined as 180° minus the angle between the two tracks in the r — (p plane. 
^The visible energy, the visible mass and the total transverse momentum of the event are calculated using 
the method described in [ p6| . 
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initial state radiation or two-photon interactions), a large overestimate of mx could occur, 
providing a high background in the accepted mass range. The following selection criteria are 
designed to reject these backgrounds. 

Bl To reject background events with two tracks and neutrinos, the acoplanarity angle of the 
two tracks is required to be ^acop < 1°- 

B2 To reduce the background from two-photon interactions and radiative events, the total 
visible energy of the event is required to be E^is > O.Q^/s and the total momentum along 
the beam direction is required to be < 0.2y/s. 

B3 The polar angle of the missing momentum vector must satisfy | cos 6miss\ < 0.8. The total 
energies measured in the forward detectors (FD, SW, GC) must satisfy: i^po < 5 GeV, 
Esw < 5 GeV, Eqq < 5 GeV. The total energy deposited in the HCAL is required to be 
less than 15 GeV. The last two requirements introduce relative inefficiencies of at most 
3.2% and less than 1%, respectively, due to electronic noise in these detectors; the final 
efficiencies are corrected by these factors. 

B4 Both selected tracks are required to be at least 0.5° away from the CJ anode planes in 
the r — (f) projection, to avoid momentum mismeasurements. 

At this stage the background is mainly due to /i^/x^ events with a smaller contribution from 
r~^r~ events where both r's decay into muons. The intermediate-mass signal is then selected 
by requiring that: 

B5 For both tracks, m| > {0.27y/sy. At v^=130-136, 161 GeV this cut is not stringent 
enough, therefore it is required that m'^ > (45 GeV)^. Figure ^(b) shows the distribution 
of mx for both selected tracks, after cut B4 at ^/s = 183 GeV. 

The effect of cuts B1-B5 for data and simulated events can be seen in Table |l|. 

An event is selected as a candidate if it satisfies either of the two sets of selection criteria 
(A1-A2) or (B1-B5), following the preselections (P1-P5). The overall detection efficiency for 
spin-0 particles is shown in Figure 0(a), as a function of mx, for y/s = 183 GeV. It exceeds 
37% even in the mass interval where the dE'/dx-based selection is inefficient. At lower centre- 
of-mass energies the efficiencies have similar values and behaviours. For spin-1/2 particles, the 
efficiencies are 2-9% lower due to the different angular distribution of the tracks. 

After this selection one candidate survives in the ^/s = 161 GeV data set while no candidates 
are found in the other data samples. The expected backgrounds at 130-136, 161, 172, 183 GeV 
energies are 0.46, 0.22, 0.31 and 1.01 events, respectively (see Table |]). 

The candidate at ^/s = 161 GeV has selected tracks with momenta pi = (51.8 ± 3.8) GeV 
andp2 = (46.7±2.9) GeV, leading to estimated masses of (61.6±3.2) GeV and (65.6±2.1) GeV, 
respectively. A possible interpretation for this event in the SM is a t~^t~ event where both r's 
decay into /x's. 
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Cuts 


Data 


Background Simulation 


Signal MC (%) 






Total 


e"'"e 








Others 




^55 


^80 


Pl-2 


1595 


1517.04 


806.97 


378.71 


90.58 


147.05 


93.73 


97.0 


97.5 


97.7 


P3-4 


1186 


1158.97 


636.45 


309.29 


66.48 


109.11 


37.63 


92.7 


93.1 


95.8 


P5 


369 


357.11 


0.00 


288.82 


5.07 


56.53 


6.70 


92.5 


92.8 


95.4 


A1,A2 





0.00 


0.00 


0.00 


0.00 


0.00 


0.00 


84.0 


7.5 


95.4 


Bl 


297 


284.40 


0.00 


247.73 


2.55 


32.30 


1.81 


88.8 


90.2 


93.0 


B2-4 


88 


92.61 


0.00 


92.30 


0.25 


0.00 


0.06 


61.1 


61.4 


0.2 


B5 





1.01 


0.00 


0.85 


0.15 


0.00 


0.01 


8.7 


32.0 


0.0 


Comb. 





1.01 


0.00 


0.85 


0.15 


0.00 


0.01 


84.1 


37.5 


95.4 



Table 1: The numbers of events remaining after each cut for data collected at ^/s = 183 GeV 
and for various Monte Carlo background processes normalised to the integrated 
luminosity of the data ("Others" refers to e"'"e~ — > qq, and e~^e~ four-fermion 
processes) . In the last three columns the efficiencies for are given (in percent) 
for mx = 45, 55, 80 GeV. The last line gives the results after combining the two 
analyses described in the text. The relative rate of data and backgrounds and the 
composition of the background at lower energies is similar to this case. 



4.2 Systematic uncertainties 

The main systematic errors affecting the signal detection efficiencies are listed below and their 
values are summarised in Table |^. These are: the statistical error from the Monte Carlo 
samples, the uncertainty due to the linear interpolation of the efficiencies, the evaluation of the 
integrated luminosity |2^ , the measurement of 0acop and the measurement of E^is and . The 



uncertainty due to the dE/dx measurement (cuts Al and A2) is estimated as follows: for each 
Monte Carlo signal event, if both tracks have dE/dx > 11 keV/cm [dE/dx < 9 keV/cm), both 
dE/dx values are decreased (increased) by their errors, then the selection cuts are applied to 
the modified event. The resulting efficiency is therefore lower; the variation in the efficiency is 
taken as the systematic error. The uncertainty due to the momentum measurement (cut B5) is 
estimated in the same way, but instead of varying the dE/dx values, both track momenta are 
increased by their errors. Again, the resulting decrease in efficiency is taken as the systematic 
error. The uncertainty due to cut B5 is estimated to be at most 30% in the intermediate mass 
region, where it dominates all the systematic effects. 

The systematic errors are assumed to be independent, and the total systematic error is 
calculated as the quadratic sum of the individual errors. 

The uncertainty on the background arises essentially from the cuts Al, A2 and B5 and 
is studied in the same way as described above. At ^/s = 183 GeV the uncertainty on the 
background is ± 0.69 events. In computing the results, no background subtraction is performed. 



4.3 Search for particles with charge ±2/3 

The search for stable and long-lived massive particles with charge ±2/3 has been performed 
without assumptions for the interaction properties of such particles. Therefore, requirements 
on energy deposits in the calorimeters are not applied. The selection is based on dE/dx 
information only. The applied cuts are PI through P4, followed by Al and A2. 
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Quantity 


Systematic uncertainty (%) 




{mx/Vs < 0.27 or mx/V^ > 0.36) 


(0.27 < mx/Vs < 0.36) 


MC statistics 


0.9-1.6 


Luminosity 


0.5-0.9 


0acop 


0.1 


E ■ 


0.5-0.9 


Interpolation 


2.0 


5.0 


dE/dx 


0.0-1.8 


1.8-8.9 


mx 


0.0-1.7 


1.7-30.0 


Total 


2.1-3.5 


5.7-30.0 



Table 2: Relative systematic uncertainties associated with the various quantities used. The 
systematic errors vary shghtly with centre-of-mass energy, but strongly with mx- 
The uncertainty in the nix determination dominates in the intermediate mass region 
but has a small effect in the mass regions where the dE/dx-based selection is efficient. 



The ionization energy loss dE/dx scales roughly as (Q/e)^; hence particles having \Q/e\ = 2/3 
would generally have dE/dx values smaller than those of charge ±1 particle tracks. For 
masses mx < O.ily/s (e.g., at ^/s = 183 GeV, mx = 75 GeV, one has p = 79 GeV and 
dE/dx = 8 keV/cm) or mx > O.U,/s (e.g., at = 183 GeV, mx = 82 GeV, one has 
p = 61 GeV and dE/dx = 13 keV/cm), the dE/dx for the two tracks of the signal events 
would be well separated from that of SM di-lepton events, as can be seen from Figure 0. The 
narrow intermediate region is not covered by any kinematical analysis. However, the efficiency 
of selections A1-A2 remains larger than 20%, because the small variation of track momenta 
due to initial state radiation can yield a variation of dE/dx which is larger than the excluded 
range. 

After this selection, one candidate remains in the data sample at ^/s = 183 GeV, and no 
candidates are left in any of the other data sets. The background at i/i = 183 GeV is estimated 
to be less than 0.3 events; the total background, summed over all energies, is estimated to be 
less than 0.45 events, as shown in Table |^. The two tracks of the candidate event have dE/dx 
values of 11.06 ± 0.35 keV/cm and 11.06 it 0.34 keV/cm, and momenta (70.9 ± 7.0) GeV and 
(64.7±6.6) GeV, yielding estimated masses of (93.4±9.6) GeV and (88.0ib9.8) GeV, respectively 
(the particle mass is evaluated from the momentum and the estimate of jSj obtained from the 
dE/dx measurement, assuming \Q/e\=2/3). However in this event the mismatch in polar angle 
between the two selected tracks and the ECAL clusters gives evidence of a reconstruction error 
of both tracks. This effect is described by the OPAL detector simulation program and would 
explain an overestimate of the dE/dx values by affecting the path length. 

The systematic uncertainties on the selection efficiencies are estimated as described in the 
previous section and they are similar in size. 

5 Results 

The numbers of candidates found in the search for particles with charge |Q/e| = l and \Q/e\=2/3 
are summarised in Table ^, together with the expected backgrounds. The data show no excess 
above the expected background from Standard Model processes. Therefore model-independent 
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cross-section upper limits have been computed for the pair-production of massive charged long- 
hved particles. 









\Q/e\ = 


1 search 






= 2/3 search 




Lumi. 


candidates 


background 


cand. 


back. 


(GeV) 


(pb-i) 


dE/dx kin. 


comb. 


dE/dx kin. 


comb. 






130-136 


11.0 








0.03 0.42 


0.46 





<0.08 


161 


10.1 


1 


1 


<0.005 0.22 


0.22 





<0.02 


172 


10.3 








<0.005 0.31 


0.31 





<0.07 


183 


58.1 








<0.005 1.01 


1.01 


1 


<0.27 


Total 


89.5 


1 


1 


0.03 1.96 


2.00 


1 


<0.44 



Table 3: The number of candidate events and the expected background at all energies, for 
the search for \Q/e\ = 1 and 2/3 particles. For \Q/e\ = 1, the result of both the 
dE/dx-based and the kinematic selections are shown separately and then combined. 
In the second column, the integrated luminosity is given for each energy. 



All data collected so far at energies above the Z° peak {y/s = 130 — 183 GeV) have been 
combined, assuming s-channel production and therefore an energy dependence of the cross- 
section of for spin-0 particles and |(1 — ^) for spin-1/2 particles, where /3 = p/E ~ 

1 — Am'j^/s. 



For the production of particles of charge ±1, the candidate event at = 161 GeV has 
been taken into account, assuming a mass of (64.4 ± 1.8) GeV, equal to the weighted average of 
the two measured masses, mentioned in Section |4.1| . Similarly, for particles of charge ±2/3, the 
candidate event at y/s = 183 GeV is included with a mass of (90.7 ± 6.9) GeV. In evaluating 
upper limits, the candidates are counted in mass intervals centred on their central values and 
±2(7 wide. The total systematic error is incorporated into the limits, following the prescription 



of Ref . p8| . No background subtraction is performed. 

In Figure |^, the 95% CL upper limit on the cross-section at y/s = 183 GeV is shown for 
spin-0 particles. The 95% CL upper limit on the pair-production cross-section varies from 0.05 
to 0.1 pb in the mass range 45 < mx < 89.5 GeV. The cross-section limits are compared with 
the predicted cross-sections for pair-production of right- and left-handed smuons and staus. 
For these two slepton species, the production cross-section does not depend on the MSSM 
parameters but only on the slepton mass. The 95% CL lower limits on the mass of right- and 
left-handed smuons and staus of 82.5 GeV and 83.5 GeV, respectively, are derived, as shown 
in Figure |^. 

Figure §(a) shows the 95% CL upper limit on the cross-section at a/s = 183 GeV for spin-1/2 
particles; the limit varies from 0.05 to 0.19 pb in the mass range 45 < mx < 89.5 GeV. This limit 
is compared with the predicted cross-sections for chargino production and heavy charged lepton 
production. The MSSM parameters have been chosen to minimise the predicted chargino cross- 
section at every chargino mass value (assuming a heavy sneutrino, mj> > 500 GeV), without 
any restriction on the mass of the lightest neutralino. Therefore a 95% CL lower limit on the 
masses of long-lived charginos is derived at 89.5 GeV for every choice of the MSSM parameters. 
The 95% CL lower limit on the heavy charged lepton mass is 89.5 GeV. 

Figure §(b) shows the 95% CL upper limit on the cross-section at ^/s = 183 GeV for spin- 
1/2 particles of charge ±2/3; the limit varies between 0.05 to 0.2 pb in the mass range 45 < 
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mx < 89.5 GeV. For spin-0 leptoquarks the cross-section upper limits would be slightly better 
than those calculated from spin-1/2 particles. These limits can be compared with leptoquarks 
and deconfined quarks models which provide cross-section predictions. 

The results obtained are valid for particles with a lifetime longer than 10^^ s. This is 
calculated in the worst case, e.g. the heaviest (and therefore slowest) particles excluded by this 
search, requiring that the decay probability of these particles at a flight distance larger than 
3.0 m is greater than 95%. For lower mass values the results are also valid for shorter hfetimes. 

6 Summary and Conclusions 

A search has been performed for pair-production of stable and long-lived massive particles with 
charge |Q/e| = lor2/3. The primary tool used in this search is the precise dE/dx measurement 
provided by the OPAL jet chamber. No evidence for the production of heavy particles with 
masses in the range of 45 to 89.5 GeV was observed. For s-channel production, the upper 
limits on the cross-section vary between 0.05 and 0.19 pb in the range of masses explored for 
particles of charge ±1. Within the framework of the MSSM, lower mass limits on the the right- 
(left-) handed smuons and staus of 82.5 GeV (83.5 GeV) have been obtained. Heavy long-hved 
charged leptons and long-lived charginos with masses smaller than 89.5 GeV are excluded. For 
particles with charge ±2/3 the upper limits on the production cross-section vary between 0.05 
and 0.2 pb in the range of masses explored. The above limits are valid at the 95% CL for 
particles with lifetimes longer than 10~^ s. 
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Figure 1: (a) The distribution of the ionization energy loss, dE/dx, as a function of the appar- 
ent momentum p/Q for data collected at ^/s = 183 GeV. The two hatched regions 
are the search regions at -^/s = 183 GeV. The momentum lower limit is defined by 
the preselection cut p > 0.1^5. A cutoff of pt > 0.1 GeV is made to reject low 
momentum tracks curling in the jet chamber volume. 

(b ) Expanded view of the search regions. The theoretical curves for heavy long-lived 
particles are shown. In e"'"e~ X+X~ events, the momentum of the particles 
of a given mass is fixed by ^/s. For \Q/e\ = 2/3, mx = 85 GeV only the position at 
y/s = 183 GeV is visible, while at y/s =172 GeV the dE/dx value lies outside the 
plot; y/s = 161 GeV is below the pair-production threshold for this mass. 
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Figure 2: (a) Selection efEciency as a function of the mass nix for \Q/e\ = 1 and 2/3 particles, 
at y/s = 183 GeV. The dashed-dotted and dashed lines represent the efficiency of 
the dE/dx-based and kinematic selections, respectively, in the case of the spin-0 
particles of charge ±1. The solid line refers to the combined selection. The dotted 
line represents the overall efficiency of the search for spin- 1/2 particles of charge 
±2/3. 

(b) The distribution of the m^, after cut B4, for data collected at y/s = 183 GeV 
(points with error bars), for all simulated backgrounds (grey histogram) and for a 
simulated signal (dashed line) with mx=55 GeV. The signal histogram is arbitrarily 
normalised. The cut B5 is indicated by the arrow pointing to the accepted region. 
No event is selected when both selected tracks are required to pass this cut. 
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Figure 3: Model-independent 95% CL upper limits on the pair-production cross-section of 
spin-0 heavy long-lived particles of charge ±1 as a function of their mass (solid 
line). In calculating the upper limit, the candidate is considered as described in the 
text. The MSSM predicted cross-sections for right-handed (dotted line) and left- 
handed (dashed line) smuons and staus are also shown. The 95% CL lower limits 
on the masses of these sleptons are at the crossing point between the experimental 
and theoretical curves. The grey region is kinematically inaccessible. 
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Figure 4: (a) Model-independent 95% CL upper limits on the pair-production cross-section of 
spin- 1/2 heavy long-lived particles of charge ±1 as a function of their mass (solid 
line). In calculating the upper limit, the candidate is considered as described in 
the text. The MSSM predicted cross- sections for charginos (dashed line) and heavy 
leptons (dashed-dotted line) are also shown. The 95% CL lower limits on the masses 
of these particles are at the crossing point between the experimental and theoretical 
curves, (b) Model-independent 95% CL upper limits on the pair-production cross- 
section of spin- 1/2 heavy particles of charge ±2/3 as a function of their mass (solid 
line). In both plots the grey region is kinematically inaccessible. 
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